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Abstract 

We consider motion of an overdamped Brownian particle in a washboard potential exerted to a 
static tilting force. The bias yields directed net particle motion, i.e. a current. It is demonstrated 
that with an additional time-delayed feedback term the particle current can be reversed against the 
direction of the bias. The control function induces a ratchet-like effect that hinders further current 
reversals and thus the particle moves against the direction of the static bias. Furthermore, varying 
the delay time allows also to continuously depreciate and even stop the transport in the washboard 
potential. We identify and characterize the underlying mechanism which applies to current control 
in a wide temperature range. 

PACS numbers: 05.40.-a, 02.50.Ey, 05.60.-k 



1 



The dynamics of particle motion in periodic potentials at finite temperature is an exten- 
sively studied field [jj . Being exerted to a static tilt the system responds with a permanent 
transport in the direction of the bias. A number of experimental situations can be described 
on the basis of one-dimensional particle motion in a tilted spatially periodic potential 
[l^ . Here we consider particle motion in a washboard potential which is often employed 
as as a paradigm to model transport in one-dimensional periodic and symmetric structures 
-llSH. Our aim is to demonstrate that the direction of the current can be controlled 
with the application of a time-delayed feedback method. Although the delayed feedback 



method was originally proposed by Pyragas 
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to stabilize unstable states in deterministic 



systems meanwhile it has been facilitated in various other contexts [20]. Recently for the 
control of transport in Brownian motors a feedback strategy has been successfully utilized 
for two ratchet systems interacting through a unidirectional delay coupling I21!]. The effect 
of time-delayed feedback on the rectification of thermal motion of Brownian particles has 



been studied in overdamped ratchet systems 22|. 



We consider an overdamped Brownian particle moving in a tilted one-dimensional 
spatially-periodic and symmetric potential. The dynamics is governed by the following 
Langevin equation expressed in dimensionless form 

q{t) = -U'{q{t)) + F + m + f{t)- (1) 

The dot and prime denote differentiation with respect to time and coordinate respectively. 
The potential is given by 

U{q) = f/(g + 1) = - cos(27rg)/(27r) . (2) 

F is a tilting force serving for a static bias. The particle is subjected to a Gaussian distributed 
thermal, white noise C,{t) of vanishing mean {C,{t)) = 0, obeying the well-known fluctuation- 
dissipation relation {^{t)^{t')) = 2kBT6{t — t') with fc^ and T denoting the Boltzmann 
constant and temperature respectively. The last term in ([T]) denotes a continuous feedback 
term of the form 

/(t) = i^(l-tanh[g(t-r)-g(t)]) (3) 

of strength K and with delay time r. Note the restriction < f{t) < 2K. In what 
follows we set exemplarily the value of the tilting force to F = —0.75. Motion in the biased 
periodic potential is then favored to the left. The potential barrier height is modulated by 



the static tilt as AE{F) = - - \F\[tt/2 - arcsin(|F|)]}/7r . For motion to the left 
with negative velocity v there remains a potential barrier of AE{F = —0.75) = 0.032 to be 
surmounted. This has to be compared with the barrier height for the unbiased potential, 
AE{F = 0) = I/tt ~ 0.318, which we hereafter denote by AE. 

Without feedback, K = 0, the particle forced by the thermal noise escapes from a well of 
the potential into an adjacent one (due to the negative tilt preferably to the left) where it 
may dwell for a while before a further escape happens. Nevertheless the long-term behavior 
reveals a negative net motion. The expression for the corresponding average dwell time, 
t dwell, exhibits the typical exponential dependence on the ratio between the barrier height 
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and the thermal energy according to tdweii oc ex^{AE /ksT) 

Interestingly, this situation of net motion in one exclusive direction may change com- 
pletely imposing the Langevin dynamics to the time-delayed feedback. Indeed we found 
parameter constellations for which the direction of the net motion which, initially governed 
by the negative tilt, is eventually reversed due to the impact of the time-delayed feedback. 
In the numerical simulation of the system ([1]) with applied time-delayed feedback term ([3]) 
we set f{t) = in the interval t G [0,r), that is the system is affected by f{t) only for 
t ^ T. The impact of time-delayed feedback of strength K = 0.8 is studied hereafter. It is 
illustrative to introduce the effective tilt force as 

Feffit) = K{1 + tanh[g(t) - q{t - r)]) +F. (4) 

For times t < t the particle evolves in the tilted washboard potential without feedback and 
moves preferably to the left. If r is not too small (see further below) the coordinate difference 
q{t) — q{t — T) may assume in the course of time such a large negative value that the feedback 
term, as given in Eq. ([3]), is almost zero. As long as the particle continues sliding down the 
tilted washboard potential this situation cannot change. On the other hand, if the particle 
dwells for long enough time, t > r, in a potential well the coordinate difference gets not only 
constrained but may even become positive. Therefore the contribution from the feedback 
force, K{1 + tanh[g(t) — q(t — r)]), to the effective tilt force, Feff(t), exceeds sufficiently the 
remaining one coming from the static tilt force, F < 0, and motion in the positive direction 
is entailed. Provided that K > \F\ holds and if the motion in the right direction is sustained 
for a long enough period during which the coordinate difference can grow to a larger positive 
value, the bias of the washboard potential attains the value Fe//(t) = 2K + F > 0. This 
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FIG. 1: Current in dependence of the delay time with fixed feedback strength K = 0.8 and for 
different temperatures as indicated in the plot. The static tilt force is given hy F = —0.75. 

saturation level is then very likely maintained. 

Quantitatively the net particle motion is assessed by the current. The latter is measured 
by the average particle velocity which in the long-time limit is defined as 

.= limM), (5) 

t^oo t 

with (■) indicating the average over the thermal noise. 

The Langevin equation were numerically integrated by means of a two-step stochastic 
Heun solver scheme. Without feedback we found that the current obeys the relation v ~ 
—0.163 X Ethermai foi a static tilt F = —0.75 where the thermal energy is measured in units 
of the (unbiased) barrier height, viz. Ethermai = ksT/AE. 

To illustrate the features of the time-delayed feedback we show the average velocity as 
a function of the delay time r for three different fixed temperatures in Fig. [H Averages 
were performed over 500 realizations of the thermal noise. For a low thermal energy ksT = 
0.05 X AE the current virtually vanishes regardless of the delay time r as long as r > 1. 
Apparently, for delay times t > 1 the feedback restrains the difference coordinate q{t) — 
q{t — t) to small- amplitude fluctuations around zero. This in turn keeps the effective tilt 
force Feff = K + F small so that for low thermal energies barrier crossings are impeded. 
Thus with applied feedback term a transition from negative to zero current is achieved. 
This changes considerably for higher temperatures where the current can be reversed due to 
applied time-delayed feedback. More precisely for ksT > 0.1 x AE there exists an interval of 
delay times for which the current becomes positive. Moreover, as also seen in Fig. [1], within 
this interval a resonance structure arises for which it holds that the higher the thermal 
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FIG. 2: Current as a function of the feedback strength K for a fixed delay time and for different 
temperatures as indicated in the plot. The static tilt force is given by F = —0.75. 

energy the higher the maximal current and the more the position of the latter is shifted to 
a smaller Tpeak value. For each temperature there exist a critical value of the delay time 
beyond which the current becomes negative again. Eventually for large enough r values the 
current attains the level adopted for zero feedback strength, K = 0, that is the feedback 
is no longer of influence. To understand further the effect of the feedback term we plot in 
Fig. [2] the current as a function of the feedback strength K for two different values of the 
thermal energy. For fc^T = 0.2 x AE we notice that while for small but positive values of 
K the current is not modified compared to its value adopted for = there is a more or 
less sudden increase for K > 0.7 and crossover from negative to positive values of v takes 
place at K = 0.75, that is when the two contributions to the effective tilt force F^ff balance. 
Further enlargement of K leads to a comparatively rapid linear growth of the current. The 
same holds true for the lower thermal energy fc^T = 0.1 x AE except that the rapid growth 
of V at the crossover point, which is also located at K = 0.75, is not as pronounced as in 
the previous case. The current in dependence of the thermal energy Ethermai = kBT/AE, 
measured in units of the barrier height, which is shown in Fig. [3] for an exemplary fixed 
feedback strength K = 0.8 and delay time r = 16, exhibits a resonance structure. That is, 
there exist an optimal thermal energy, Etermai ~ 0.24, associated with a peak of the current 
to either side of which the current decays fairly rapidly. 

To gain insight into the role played by the delay time with regard to the current efficiency 
we represent the time evolution of the coordinate the difference q{t) — q(t — r), and the 
effective tilt force Feff(t) for a delay time r = 5 and r = 32 in Figs. H] and [5] respectively 
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FIG. 3: Current in dependence of the scaled thermal energy -Ej/iermaZ = ksT/ /S.E for fixed feedback 
strength K = 0.8 and fixed delay time r = 16. The static tilt force is given by F = —0.75. 

for one realization of thermal noise. For a delay time r ~ t dwell — 5 < Tpeafc = 32, which is 
of the order of the dwell time, the resulting current is comparatively low but positive. This 
is reflected in the temporal evolution of the coordinate in Fig. [n](a). In an early stage of 
the dynamics (prior to the application of the feedback) the particle escapes from its starting 
position (the bottom of a potential well) into the well adjacent to the left where it dwells for 
considerably long time. After another jump to the left has taken place further occasional 
jumps follow in the opposite direction. That behavior is related with the time evolution of 
the effective tilt force Feff{t) as well as the difference coordinate q{t)~q{t~T) which most of 
the time perform small amplitude fluctuations around zero apart from the occasional sudden 
bursts being correlated with the kicking of the particle from one potential well into another 
one. Nevertheless on average the particle motion proceeds to the right (cf. the long-term 
behavior of q{t) in the inset). 

In comparison for the delay time Tp^ak = 32 > tdweiu yielding the peak in the current (cf. 
Fig. [1]), the particle slides already far down the tilted potential in an early stage for t < 328. 
During this time the coordinate difference q{t) — q(t — r) stays strictly negative as seen in 
Fig.O Hence the delay term itself, K{t) = K {1 + ta.nh[q{t) — q{t — r)]), is vanishingly small 
and consequently the particle feels to the full extent the bare tilt force F < 0. Subsequently 
to a trapping interlude a longer excursion in the negative direction occurs terminated by 
yet another trapping period. For the particle dwelling inside the potential well and being 
buffeted by the thermal noise, likewise in the previous case, Feff(t) as well as q(t) — q(t — r) 
fluctuate around zero. Strikingly, in contrast to the previous case of short delay time r = 5, 
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FIG. 4: (a): Time evolution of the position for one realization of noise. The inset displays the 
long-term behavior, (b): The difference — g(t — r). (c): The effective tilt force i^e//(i)- Notice 
the assignment of the bursts in (b) and (c) to the respective particle transitions. The parameter 
values read as: F = -0.75, K = 0.8, r = 5 and ksT = 0.1 x AE. 
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FIG. 5: As in Fig. [3]for the same realization of thermal noise but for delay time r = 32. 

in the present case of r = 32 the coordinate difference retains not only a strictly positive 
but a by far larger value for a comparatively longer time which in turn serves for ongoing 
positive effective tilt force Feff(t) ^ 2K + F > 0. Conclusively, it is the interplay between 
the coordinate difference with persisting positive amplitude and the positive effective force 
that is responsible for forcing vigorously the particle in the right direction. Note that after 
current reversion the washboard potential experiences effectively a positive tilt, 2K + F > 0, 
that is larger than the absolute value, \F\, of the negative tilt started with. On the other 
hand, an increased tilt is related with a shorter average dwell time tdweii which makes it less 
likely that the coordinate difference returns to small values let alone become again negative. 
In this sense the asymmetry K > \F\ serves for a ratchet-like effect hampering a further 
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reversal of the current back from positive to negative value. 

In conclusion, we have identified and characterized a mechanism to control the thermally 
driven transport for overdamped Brownian particles evolving in a washboard potential under 
the mutual impact of static bias and time-delayed feedback. That is with appropriately 
chosen strength of the delay term and the delay time the direction of the current can be 
reversed in such a way that particle motion proceeds opposite to the starting direction 
governed by the static bias. Due to its construction the control function induces a ratchet- 
like effect that hinders further current reversals. Moreover, varying the delay time allows 
also to continuously tune the current between its possible negative and positive extreme 
value including stoppage of transport in the washboard potential. Our results apply to a 
wide temperature range. 
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